The effect of spin-orbit interaction on entanglement of two-qubit Heisenberg XYZ 

systems in an inhomogeneous magnetic field 
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The role of spin-orbit interaction on the ground state and thermal entanglement of a Heisenberg 
XYZ two-qubit system in the presence of an inhomogeneous magnetic field is investigated. For a 
certain value of spin-orbit parameter D, the ground state entanglement tends to vanish suddenly 
and when D crosses its critical value D c , the entanglement undergoes a revival. The maximum 
value of the entanglement occurs in the revival region. In the finite temperatures, there are revival 
regions in D — T plane where increasing of temperature first increases the entanglement and then 
tends to decrease it and ultimately vanishes for temperatures above a critical temperature. This 
critical temperature is an increasing function of D, thus the nonzero entanglement can exist for 
larger temperatures. In addition, the amount of entanglement in the revival region depends on the 
spin-orbit parameter. Also, entanglement teleportation via the quantum channel constructing by 
above system is investigated and influence of spin-orbit interaction on the fidelity of teleportation 
and entanglement of replica state is studied. 

PACS numbers: 03.67.Hk, 03.65.Ud, 75.10.Jm 



I. INTRODUCTION 



A. Einstein, B. Podolsky and N. Rosen, in their fa- 
mous EPR paradox, argued that in general two quantum 
system cannot be separated even if they are located far 
from each other E. Schrodinger named this quantum 
mechanical property as Entanglement 0. Today, entan- 
glement is a uniquely quantum mechanical resource that 
plays a key role in many of the most interesting applica- 
tions of quantum computation and quantum information 
[3, El- Thus a great deal of efforts have been devoted 
to study and characterize the entanglement in the recent 
years . The central task of quantum information theory 
is to characterize and quantify entanglement of a given 
system. A mixed state p of a bipartite system is said 
to be separable or classically correlated if it can be ex- 
pressed as a convex combination of uncorrelated states 
Pa and Pb of each subsystems i.e. p = Yl^iPA ® Pb 

i 

such that u>i > and Y] ujj = 1, otherwise p is entangled 

0, |. Many measures of entanglement have been intro- 
duced and analyzed 0, H, 0| , but the one most relevant 
to this work is entanglement of formation, which is in- 
tended to quantify the resources need to create a given 
entangled state :6j. For the case of a two-qubit system 
Wootters has shown that entanglement of formation can 
be obtained explicitly as: 



E( P ) = E[C(p)] = h 



(1) 



where h(x) = —xlog2X — (1 — x)log2(l 
nary entropy function and C(p) — max{0, 2A 



x) is the bi- 
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is the concurrence of the state, where AjS are positive 
square roots of the eigenvalues of the non-Hermitian ma- 
trix R = pp, and p is defined by p := {a v ®a v )p*{a v ®cr v ). 
The function 2 is a monotonically increasing function and 
ranges from to 1 as C goes from to 1, so that one can 
take the concurrence as a measure of entanglement in its 
own right. In the case that the state of the system is 
pure i.e. p = \i/j){4>\, \ip) = a\00) + 6j01) + c|10) + d\ll), 
the above formula is simplified to C(\ip}) = 2 | ad — be \. 
The spin chain is the natural candidates for the real- 
ization of entanglement and Heisenberg model is the 
simplest method for studying and investigating the be- 
havior of the spin chains . Nielsen Q is the first per- 
son who studied the entanglement of two-qubit Heisen- 
berg XXX- chain modeled with the Hamiltonian H = 
J(Ti ■ cr 2 + B • (erf + af). He showed that entanglement 
in such systems exists only for antiferromagnatic ( J > 0) 
case below a threshold temperature T c . After Nielsen, 
entanglement in the two-qubit XXX, XXZ and XY sys- 
tems in the presence of homogeneous and inhomo gen eous 
magnetic field has been investigated [jj fiol [Til. H3. fl3| . 
The effect of anisotropy due to spin coupling in the 
x,y, z direction has also studied in a number of works 
15, 16]. In ref fuj Yang et al. have shown that in XYZ 
Heisenberg systems, an inhomogeneous external mag- 
netic field can make the larger revival, improve the crit- 
ical temperature and enhance the entanglement. Spin- 
orbit (SO) interaction cause another type of anisotropy 
[H M, Hfl, [H HI [H, H. The effect of SO interaction 
in the thermal entanglement of a two-qubit XXX system 
in the absence of magnetic field has been studied in [25| . 
However, the entanglement for a XYZ Heisenberg system 
under an inhomogeneous magnetic field in the presence 
of SO interaction has not been discussed. Therefore, in 
this paper we investigate the influence of SO interaction 
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on the entanglement and entanglement teleportation of 
two-qubit system at thermal equilibrium. 
On the other hand, among the numerous concepts to 
implement a quantum bit (qubit), approaches based on 
semiconductor Quantum Dots (QDs) offer the great ad- 
vantage that ultimately a miniaturized version of quan- 
tum computer is feasible. Indeed, at first D. Loss and 
D.P. Diviencenzo proposed a quantum computer proto- 
col based on electron spin trapped in semiconductor QDs 
in 1998 [26|, H?], [28[ . Here, the qubit is represented by a 
single electron in a QD which can initialize, manipulate 
and read out by extremely sensitive devices. Compared 
with other systems such as quantum optical systems (29[ 
and nuclear magnetic resonance (NMR) [3(J EH [S3] , QDs 
are argued to be more scalable and has long decoherence 
time [3j. The above Heisenberg system is suitable for 
modeling and computing the entanglement of a two-qubit 
system represented by two electrons confined in two ver- 
tically coupled quantum dots (CQDs), respectively. Due 
to weak lateral confinement electrons can tunnel from 
one dot to another dot and spin-spin and spin-orbit in- 
teraction between the two qubits exists. In GaAs dou- 
ble QDs, hyperfine interactions have been identified to 
dominate spin mixing at small magnetic fields, while SO 
interaction is not relevant in this regime. However SO in- 
teraction and the coupling magnetic fields are expected 
to be orders of magnitude stronger in InAs compared to 
GaAs [33J. SO interaction in such nanostructures can be 
investigated with the help of quantum optical methods 

& 

Taking the advantage of tunability of SO strength [34J, 
HE l36l . [37| , we show that this type of interaction cause 
to enhancement of entanglement in the revival region, 
increase the volume of revival region and improve the 
critical values of other parameters. 

In the following, as an application of the above system, 
the entanglement teleportation of a two-qubit pure state 
via the above two-qubit system is investigated and av- 
erage fidelity between input and output states is calcu- 
lated. C. Bennet et al. have shown that two entangled 
spatially separated particles can be used for teleporta- 
tion [3g|. They also argued that states which are less 
entangled still could be used for teleportation but they 
reduce "the fidelity of teleportation and/or the range of 
state \4>) which can accurately be teleported" . After then 
S. Popscu by using hidden variable model, have shown 
that, teleportation of a quantum state via a pure classi- 
cal communication cannot performed with fidelity larger 
than I [39l | . Thus, mixed quantum channels which al- 
lows to transfer the quantum information with the fi- 
delity larger than | are worthwhile. Horodecki et al. 
have calculated the optimal fidelity of teleportation for 
bipartite state acting on C d ® C d by using the isomor- 
phism between quantum channels and a class of bipar- 
tite states and twirling operations [4(j. G. Bowen and 
S. Bose have shown that "standard teleportation with an 
arbitrary mixed state resource is equivalent to general 
depolarizing channel with the probabilities given by the 



maximally entangled component of the resource. This en- 
ables the usage of any quantum channel as a generalized 
depolarizing channel without additional twirling opera- 
tion" [4ll |. Using the property of linearity of teleporta- 
tion process [38(, J. Lee and M.S. Kim have shown that 
" quantum teleportation preserves the nature of quantum 
correlation in the unknown entangled state if the chan- 
nel is quantum mechanically correlated" and then they 
investigate ent ang lement teleportation via two copies of 
werner states [42[. The entanglement teleportation via 
thermally entangled state of a two-qubit Heisenberg XX- 
chain and XY-chain has been studied by Ye Yeo et al. 
[43I l44j . The effect of spin-orbit interaction on entangle- 
ment teleportation on a two-qubit XXX-Heisenberg chain 
in the absence of magnetic field is reported by G. Zhang 

In this paper we investigate the ability of the above men- 
tioned two-qubit system for the entanglement teleporta- 
tion. We show that spin-orbit interaction and inhomoge- 
neous magnetic field have effective effect on the entangle- 
ment of replica state and the fidelity of teleportation. A 
minimal entanglement of the thermal state in the model 
is required to realize efficient entanglement teleportation 
and we can attain to this minimal entanglement, in the 
case of J z < 0, by introducing SO interaction. 
The paper is organized as follows. In sec. 2 we introduce 
the Hamiltonian of a two qubit Heisenberg system under 
inhomogeneous magnetic field with SO interaction and 
write the thermal density matrix of the system related 
to this Hamiltonian and ultimately calculate the thermal 
concurrence of the system. In subsection 2.1, the ground 
state entanglement is calculated and results are plotted in 
figs. l-2c. In subsection 2.2 the finite temperature entan- 
glement of system is computed. The figs. 3-5 illustrate 
the obtained results. The entanglement teleportation of 
a two-qubit pure state and its fidelity derived in sec. 3 
and results are plotted in figs. 6-8. In sec. 4 a discussion 
concludes the paper. 

II. MODEL AND HAMILTONIAN 

The Hamiltonian of a two-qubit anisotropic Heisenberg 
XYZ-model in the presence of inhomogeneous magnetic 
field and spin-orbit interaction is: 

H = \{ J X of o-f + J v a\a\ + J z o\a 7 2 ' + B x • <r a + B 2 • <x 2 

+ D (<7! X(T 2 )+5 (T 1 -T-(T 2 ), (2) 

where crj = (<rj , er J , er|) is the vector of Pauli ma- 
trices, Bj (j = 1,2) is the magnetic field on site j, 
Jfj_ (/j = x, y, z) are the real coupling coefficients (the 
chain is anti-ferromagnetic (AFM) for > and ferro- 
magnetic (FM) for < 0) and D is called Dzyaloshinski- 
Moriya vector, which is first order in spin-orbit coupling 
and is proportional to the coupling coefficients (J M ) and 
r is symmetric tensor which is second order in spin- 
orbit coupling [l8|, Ojl, HH ■ For simplicity we assume 
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Bj = Bj z such that B\ = B + b and B2 = B — b, where b 
indicates the amount of inhomogeneity of magnetic field. 
The vector D and the parameter S are dimensionless, in 
system like coupled GaAs quantum dots |D| is of order 
of a few percent, while the order of last term is 10~ 4 and 
is negligible. If D = J Z D z and ignore the second order 
spin-orbit coupling, then the above Hamiltonian can be 
expressed as: 



In the above equations 



and M ± 



(b±Q 2 



are the normalization constants. Here we 



define, f := yjb 2 + J 2 + (J Z D) 2 and r] := V ' B 2 + (J7) 2 , 
for later convenience. 

The state (density matrix) of a system in equilibrium at 
temperature T is p = Z~ x exp(— r^f)j where Z is the 
partition function of the system and fee is the Boltzman 
constant. For simplicity we take fcs = 1. In the standard 
H = •^7( cr i a i + °2 ) + (J + iJzD)a l a 2 + ( J — iJzDjo^^g^ the density matrix of the system in the thermal 
J z _ , .B + b. , .B — b. , equilibrium can be written as: 

+ y^o-a + C— J- K + (~ y- K> W 



where J 



is the mean coupling coefficient 
in the XY-plane, 7 := '! x , J T y , specifies the amount 
of anisotropy in the XY-plane (partial anisotropy, 
— 1 < 7 < 1) and cr = ^(a 1 ± ct 2 ') are lowering and 
raising operators. 

The Hamiltonian ([3]), in the standard basis 
{|00), |01), 1 10) , |11)}, has the following matrix form: 



Pt = 



where 



/!+ v 

io a z 

z* w 2 

/i_ 



(8) 



(\J Z + B 



V ^7 





.1 7 





b J + iJ z D 
■J — iJ z D -\J z -b 


The spectum of H is easily obtained as 

H\^ ± ) = e l ^ ± ), 



H\S k )=e i .\S k ) 








I / 



(4) 



(5) 



where the eigenstates and the corresponding eigenval- 
ues are, respectively 



I^ ± ) = ^ ± (-(7to)I° 1 ) + I 10 ))' 

|E ± )^M ± (-(^)|00) + |11)), 



(6) 



and 



£3,4 = \ Jz ±V- 



(7) 



W12 



' 2 



(cosh /?77 + — sinh /3?y) , 



^-(cosh/3^T |sinh/3£), 



_ 0J* 

J7 e 2 



(9) 



sinh /3t7 , 

sinh (3£ . 



(J+iJ z D)e 2 



and Z = 2e^ L (cosh/3£ + e~@ Jz cosh j3rj) . 
In what follows, our purpose is to quantify the amount of 
entanglement of the above two-qubit system versus the 
parameters of the system, with the main concerning on 
D . For density matrix in the form ([8]), one can show, 
by straight forward calculations, that the square roots of 
the eigenvalues of matrix R — pp are: 



Al,2 — 



I y/WiW 2 ± I z 1 1 

Mi 
e 2 



i 2 + J 2 + (J Z D) 2 sinh 2 P£ ± yj J 2 + {J Z D) 2 sinh p£ 



(10) 



^3,4 = I v^+M- T ^ 

e 2 



rj 2 + ( J7) 2 sinh 2 /3ry + J7 sinh /3?7 
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loss of generality we can assume J > and 7 > 0, since 
the above formula are invariant under substitution J — > 
—J and 7 — > —7. For the special case D = 0, these 
equations give the same results as ref [TtJ • The obtained 
results are given in the following subsections. 



Ground state entanglement and critical 
parameters 



FIG. 1: The ground state concurrence vs. b for D = (solid The behavior of the system at the quantum phase tran- 

line), D = 0.8 (dashed line) and D = 1 (doted line) where sition (QPT) point [49|], such as B = B c , b = b c or 
B — 0.8, J = 1, J z = — 1 and 7 = 0.5. D = D c , may be determined from the density matrix 

at T = 0, at which the system is in its ground state. If 
we consider T — > (or (3 — > 00) then the concurrence C 
Now, it is easy to calculate the concurrence. Without can be written analytically as: 



C(T = 0) = < 



r 



if £,<V~Jz 

if £ = rj- J z 
if i>r]-J z 



(11) 



Also, this formula can be derived by calculating en- 
tanglement of ground state of H directly: When £4 < e 2 
(or £,<!]— J z ) the ground state of H is and then 

C = I ^j- I and when e 2 < £4 (or £ > 77 — J z ) the ground 

state of H is |V>~) and then C =| V J2 +^ g ) 2 |_ 0n the 
other hand, at critical point (where £4 = £2 or £ < r]—J z ), 
the ground state of the system is an equally mixture of 
|XT) and |V>~) i.e. |G5) = (|XT) + |^ _ ))/2 and then 



C = 



The concurrence C as a function of & (at T = ) for three 
values of D (D = 0,0.8,1 ) are plotted in fig. 1. With 
increasing 6, the concurrence C is initially constant and 
equal to C =| ^ |= 0.53, then drops suddenly at critical 

value of b (& c = v /(^-^) 2 -((^£>) 2 + ^ 2 ))- At this 
point (T = 0, b = b c ), the concurrence becomes a non- 
analytical function of 6 and QPT occurs [13 • For & > b c , 
concurrence C undergoes a revival before decreasing to 
zero. The amount of concurrence at revival region de- 
pends on D, by increasing D the revival is greater. Fur- 
thermore, by increasing D, b c decreases i.e. for larger D, 
the critical point and hence revival phenomenon occurs 
in smaller b. The role of D is more obvious in figs. 2a-2c 
where concurrence is plotted in terms of D. These figures 
show that, when D reaches its critical value defined by 



Dr. 



y/(v~ J Z ) 2 -(V + J 2 ))/\U 



(12) 



the concurrence drops and exhibit a revival phenomenon 
when D crosses its critical value D c . In the revival region 
(larger D) the concurrence reaches its maximum value 
(C — 1). Fig. 2a illustrates the ground state concurrence 
variation versus D for three values of b, by decreasing b, 
the critical value of D {D c ) increases. Fig. 2b shows 
ground state concurrence versus D for three values of 7, 
by increasing 7, D c increases (amount of 7 determines the 
value of entanglement before reaching the critical point). 
In fig. 2c the ground state concurrence versus D is plot- 
ted for three values of J 2 , when | J z | increases, D c also 
increases. 

Hence, we can control value of D c by adjusting the pa- 
rameters of the system such as 6,7 and J z . We demand 
to decrease D c (because for D > D c concurrence will be 
maximized). Therefore we should choose, | J z \= 1,7 
small and b as large as possible. 



B. Thermal entanglement 

Since the relative magnitude of A,(z = 1,2,3,4) 
depends on the parameters involved, they cannot be 
ordered by magnitude without knowing the values of 
the parameters. This prevents one from writing an 
analytical expression for concurrence. For particular 
parameters, C can be evaluated numerically. The 
role of each parameter can be seen by fixing other 
parameters and drawing variation of entanglement for 
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FIG. 2: The concurrence of the ground state vs. D. a) for 
6 = (solid line), b = 0.5 (dashed line) and b = 1 (doted 
line), J z = -1 and 7 = 0.8 b) for 7 = 0.9 (solid line), 7 = 0.5 
(dashed line) and 7 = (doted line), J z = — 1 and b = 0.75 c) 
for J z — — 1 (solid line), J z = 0.5 (dashed line) and J z = 0.1 
(doted line), b = 0.75 and 7 = 0.7. 



specific values of that parameter. Cross influence of two 
parameters can be shown in 3D plots of entanglement. 
Thermal entanglement versus system's parameters is 
depicted in figures 3-5. The concurrence as a func- 
tion of T and D is shown in fig. 3. The analysis of 

results of this figure are given in the following. Let 

4 

C'j = 2Xj - (j = 1,2,3,4), it is evident that the 

i=l 

function Cj = max{0,C'j} is the concurrence of the 
system if and only if Xj — A max = max{Ai,A4} (Since 
J, 7 > we have A4 > A3 and Ai > A2). We can divide 
the regions of fig. 3 in four parts: 

i) Region for which T < T c \(see below) and D < D c \ 
in this region X max = A4 and then C4 = A4 — A3 — A2 — Ai 
determines the amount of entanglement. According to 
equation (|10p. the values of A4 depends on 7 and rj (or 
equivalently 7 and B ). Thus in this region we can ma- 
nipulate the amount of entanglement by adjusting 7 and 



FIG. 4: Thermal concurrence vs. T for different values of D. 
The parameters are the same as fig. [3] In this case D c ~ 4.51. 



B. The first critical temperature (T c i ), is the point at 
which A max = A4 = Ai and hence C = 0. 
ii) Revival region, for which T c \ < T < T C 2 (see below) 
and D < D c ; in this region X max = X\ and hence C[ 
determines the amount of entanglement. According to 
equation (flU|) . the value of Ai is adjustable by changing 
values of D and £ (or equivalently D and b) and hence 
the parameters D and b play an important role in quan- 
tifying the amount of entanglement, in this region. For 
T > T c i, the value of Ai and also the rate of enhance- 
ment of the function Ai with T increases as D increases. 
Enhancement of the function Ai with T, cause C' x rise to 
a positive number and thus the entanglement undergoes 
a revival. Since the rate of enhancement of the function 
Ai with T is an increasing function of D, the amount 
of revival increases as D increases. When T reaches the 
value T = T C 2 (second critical temperature), Ai tend to 
zero again and thus the entanglement vanishes. 
Hi) Region for which D > D c , for all values ofT; in this 
region X max = Ai and C[ is the entanglement indicator. 
The maximum value of the entanglement occurs in this 
region. At zero temperature the entanglement has its 
maximum value and by increasing the temperature the 
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FIG. 5: (color online) (a) Thermal concurrence vs. b and D, 
where B=5, (b) Thermal concurrence vs. B and D, where 
b=2. Here T = 0.3, J = 1,J Z = 0.5 and 7 = 0.3 



system loses its entanglement and ultimately vanishes at 
T = T c2 . In this region no revival phenomenon occurs. 
iv) Region for which T > T c2 for all values of D; in this 
region all values of C'- (j — 1,2,3, 4) have negative values 
and then the entanglement is zero for all values of D and 
the other parameters. 

Notice that, T c \ and T c2 are sensitive functions of D. In 
fig. 4, we try to demonstrate these facts, by illustrating 
the function of thermal concurrence vs. T for few values 
of D. This figure shows that for D < D c , T c \ is a de- 
creasing function of D but for D > D c , T c \ is undefined. 
In contrast, T c2 is an increasing function of D for all val- 
ues of T, i.e. we can create and maintain the nonzero 
entanglement at larger temperatures. In summary, fig. 
3 shows that: For T cl < T < T c2 and D < D c , there 
are regions in. D — T plane where increasing of tempera- 
ture first increases the entanglement (revival region) and 
then tends to decrease the entanglement and ultimately 
for T > T c2 entanglement vanishes. The maximum en- 
tanglement exists at zero temperature and for large D. 
In the revival region and region for which D > D c , in- 
creasing D causes the entanglement to increase. In this 
region, T c \ decreases and T c2 increases as D increases 
and hence the width of the revival region increases. In 
all regions T c2 is an increasing function of D, thus when 
D is large enough, the entanglement can exists for larger 
temperatures. Furthermore, the parameter D plays the 
role of parameter b. Fig. 5a, shows the variation of en- 
tanglement as a function of b (inhomogeneity of magnetic 
field) and D. For fixed D, there are three region in this 
figure i) Main region where b < b c ; in this region entan- 
glement is constant ii) Collapse region where b = b c ; in 
this region entanglement decreases suddenly Hi) Revival 
region where b > b c \ in this region entanglement under- 
goes a revival. Amount of D determines b c and hence 
edge of revival region. Furthermore, in the revival re- 
gion increasing D increases the entanglement. In fig. 5b, 
thermal entanglement is plotted vs. magnetic field (B) 
and D. The role of D is similar to its role in fig. 5a, 
enhancement of D improves B c and increases amount of 
entanglement in the revival region. 

One can use the above entangled two-qubit system for 
performance the teleportation protocols. The next sec- 
tion is spend to this subject. 



III. THERMAL ENTANGLEMENT 
TELEPORTATION 

For the entanglement teleportation of a whole two- 
qubit system, a thermal mixed state in Heisenberg spin 
chain can be considered as a general depolarizing chan- 
nel. Now we consider Lee and Kim's two qubit telepor- 
tation protocol (Pi), and use two copies of the above 
two-qubit thermal state, pr <8> pr, as resource Simi- 
lar to standard teleportation, entanglement teleportation 
for the mixed channel of an input entangled state is de- 
stroyed and its replica state appears at the remote place 
after applying a local measurement in the form of linear 
operators. We consider as input a two-qubit in the spe- 
cial pure state \ip m ) = cos0/2|lO) + e^sin0/2|Ol) (0 < 
< ir,0 < cf> < 27r). The density matrix related to \ipi n ) 
is in the form: 



/o 

a c 
c* b 
\ 



where a — sin 2 6/2, b = cos 2 0/2 and c = 



(13) 



-i(f> 



sm( 



Therefore concurrence of initial state is Cj„ = 2 | 
e 1 ^ sin 8/2 cos 9/2 \= sin9. The output (replica) state 
Pout can be obtained by applying joint measurement and 
local unitary transformation on input state pi n . Thus the 
out put state is given by [4l| 



Pout 



\i,V 



(14) 



where p,v = 0,x,y,z (a = I), p M „ = 
tr[E^ p c hannei\tr[E' J p channel] such that J2 V = 1 and 

u,u 

Pchannei represent the state of channel which used for 
teleportation. Here E° = |*-)(*-|, E 1 = |$-)($-|, 
E° = |$+)($+ 1 and E° = |*+)(*+| where 1^) = 
and |$±) = (|oo)±|ii)) 

arc the Bell states. 
By considering the two-qubit spin system as a quantum 
channel, the state of channel is p c hannei — Pt given in 
the equation Q and hence one can obtain p ou t as 




where 



a 



Pout 



(u>i + w 2 ){p + + p ), 



(15) 



k = 4 Re[z] v cos (j) sin 0, 
a 1 = (p + + p~) 2 cos 2 | + (wi + w 2 ) 2 sin 



2 



(16) 



b' = (wi + w 2 ) 2 cos 2 | + (p+ + p~) 2 sin 2 



d = 2e"^((Re[z]) 2 + e 2l ^v 2 ) sin 
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Following the above we can determine concurrence of out PoutPout, he. A.^s. It is easy to show 
put state by calculating positive square roots of R ou t = 



C 2 „(cosh 2 /3t7 - e 2 ?- 1 * cosh 2 f3£) 2 + 4e 2 /"* cosh 2 /3t? cosh 2 /3£ 



± 



2(cosh/3?7 + C 3 - 7 - cosh /?£) 2 
|C m ((4l) 2 e 2 ^sinh 2 /3r;+ (f) 2 e 2 ^ sinh 2 /^)| 



^3,4 — 



2(cosh/3?7 + cosh/3^) 2 

'2' 



cosh /3fy cosh /3£ ± Ci„ ( -2 ) ( j ) sinh /3?7 sinh /3£ cos < 



(cosh/3?7 + e^ Ja cosh 2 
I 



(17) 
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Thus C ouf = C(p ou t) = max{0, 2A' max - J2 K} is 

i=l 

computable when the parameters of channel are known. 
The function C ou t is dependent on the entanglement 
of initial state and the parameters of the channel 
(which determine the entanglement of channel). The 
C ou t is nonzero only for particular choice of channel's 
parameters for which Channel is greater than a critical 
value. The figs. 6 and 8 depict behavior of C ou t versus 
the parameters of the channel and Ci n . For the case 
J z > 0, the entanglement of replica state (C out ) increases 
linearly as Ci n increases. The rate of this enhancement 
is determined by D (indeed, by C 'channel (D)). But for 
the case J z < 0, C ou t is zero for small values of D. 
As D crosses a threshold value, C ou t increase when 
Ci n increase with the rate determined by amount of D 
(equivalently C 'channel (D)). The figs. 6b and 6c show 
that the parameter of inhomogeneity (b) can plays role 
ofD. 

-The Fidelity of entanglement teleportation: Fidelity 
between pi n and p out characterizes the quality of tele- 
ported state Pout- When the input state is a pure state, 
we can apply the concept of fidelity as a useful indica- 
tor of teleportation performance of a quantum channel 
1481 ] . The maximum fidelity of pi n and p ou t is defined 
to be 




F(pin,Pout) = {Tr[V {ptn) 2 Pout(Pin) 2 }} 2 

= (lpm\Pout\lpin)- (18) 

By substituting pi n and p ou t from above, we have 

F(p m ,Pout) = a' sin 2 f + b' cos 2 f + Refc'e^] sin (9(19) 

simplifying the above formula we find that the 
maximum fidelity F(pi n ,p ut) depends on initial 
entanglement (Ci n ) : 



F(p m ,Pout)=f c +f q C? n 



(20) 



FIG. 6: (color online) Entanglement of output state (Cout) 
vs. the channel's parameters and d n . Cout vs. D and d„ 
for a) 6 = 0.5 and J z > and b) b = 0.5 and J z < 0. Also, 
Cout vs. b and D for c) Cm = 1 and J z > and d)d n = 1 
and J z < 0. Where J = 1, T = 0.1, B = 1 and 7 = 0.3. 



where f c = (wi + W2) 2 and f q — | — (w\ + w 2 ) + 
2(y 2 cos 20 + (Re[z\) 2 ). The functions f c and f q are de- 
pendent on the channel's parameters only (we consider 
<f> = 0) and hence are relate to the entanglement of the 
channel. This formula is the same as the results of ref. 
42] [50], but in spite of the werner states, j q can be a 
positive number for Heisenberg chains. This means that, 
there exists a channel such that it can teleport more en- 
tangled initial state with more fidelity! but it is not a 
useful claim because, when we choose the parameters of 
the channel such that f q >0 then f c decreases and ulti- 
mately F(pi n , pout) becomes smaller than | which means 
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FIG. 7: (color online) F(p ir> 
and b) J z < 0. Where J 
7 = 0.3 



p ou t) vs. C in and D for a) J z > 
= 1,T = 0.1, -B = 1,6 = 0.5 and 



that the entanglement teleportation of mixed state is in- 
ferior to classical communication. Thus, to obtain the 
same proper fidelity, the larger entangled channel are re- 
quired for larger entangled initial state. The fig. 7 em- 
phasize the above notes. 

The average fidelity Fa is another useful concept for char- 
acterizing the quality of teleportation. The Average fi- 
delity Fa of teleportation can be obtained by averaging 
F(pi n , Pout) over all possible initial states: 



F A 



J 27r d(j> F(p out , p in ) sin 9d9 



4tt 



(21) 



and for our model Fa can be written as 



F A 



cosli 0T) + e 2f)J * (2£ 2 cosli /3£ + J 2 sinli (3£) 
3£ 2 (cosh Jh) + eP J > cosh/3^) 2 



(22) 



In the case of an isotropic XXX Heisenberg chain in the 
absence of magnetic field with spin-orbit interaction, this 
equation gives the same as results of ref. [45j • The func- 
tion Fa is dependent on the channel's parameters. The 
fig. 8 gives a plot of Fa, C out , F(p in ,p out ) and Cchannei 
versus the channel's parameters . This figure shows that, 
in the case of J z > 0, Fa, F(pi n , pout) and C ou t decrease 
when D increases. In this case Fa approaches | for large 
values of D. 

In contrast, in the case of J z < and for small values 
of Z?, Fa and F(pi n , p out ) has a constant value (smaller 
than |) and C ou t is zero. As D becomes larger than a 
threshold value, Channel undergoes a revival and then 
decreases for larger values of D. Since in the revival re- 
gion, C c hannei has its maximum value, F A , F(p mi p out ) 
and C ou t increase in this region such that for particular 
interval of D, Fa becomes larger than | and ultimately, 
tends to | for larger D. In summary, the fidelity of tele- 
portation and entanglement of replica state are depen- 
dent on the entanglement of channel which is tunable by 
the channel's parameters (such as D,b, B, ...). The ef- 
fect of D is more desirable in the case of J z < 0. In 
this case for a certain values of D, Fa becomes greater 
than | , this make the channel useful for performance the 
teleportation protocol. For large D, Fa tends to |. 




FIG. 8: (color online) a) Average fidelity (a), Output entan- 
glement (b), Fidelity between pt n and p ou t (c) and the entan- 
glement (d) of channel versus J z and D. Where J = 1,T = 
0.1, B = 1, b = 0.5, 7 = 0.3 and C in = 1 



IV. DISCUSSION 

The entanglement of a two-qubit XYZ Heisenberg sys- 
tem in the presence of an inhomogeneous magnetic field 
and spin-orbit interaction is investigated. By turning on 
the spin-orbit interaction we can change the behavior of 
the system without manipulating the other parameters. 
We have shown that the critical values of B, b and T are 
adjustable by D and we can improve the critical values of 
B, b and T. Increasing D cause to volume enhancement 
of revival region and also enhancement of entanglement 
in the revival region. Also, entanglement teleportation 
via two copies of above two-qubit system is studied. We 
have shown that, by introducing SO interaction, the en- 
tanglement of replica state and fidelity can be increased 
for the case of J z < 0, in spite to the case J z > 0. When 
D becomes very large, the fidelity approaches |, which 
is the maximal value for the classical communication. 
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